h//)*SS?C£- / 33^3 77 


nasa-cr - 173361 

^9840014932 


A Reproduced Copy 

OF 

/ifftjft 377.-/73,37? 


l-ANGLEY RESEARCH CENTER 
LIBRARY, NASA 
HAMPTON,. VIRGINIA 


Reproduced for NASA 
by the 

NASA Scientific and Technical Information Facility 



FFNo 672 Aug 65 



A>«/f 



■ / ,c: 
/ ^ 


(E34-10098) SPECTEORADICMETRIC CALIBRATION 
OR THE THEMATIC MAPPER AN E MULT ISPECTR AL 
SCANNER SYSTEM Quarterly Report, 1 Nov. 

1 9d3 - 1 Feb. 1984 (Arizona Univ., Tucson.) 

35 p HC A 03 /Mi' A01 CSCI 0BB G3/U3 


N84-2 30 00 
THRU 

N84-23002 

Unclas 

00098 



FIFTH QUARTERLY REPORT 
ON 


"SPECTRORADIQMETRIC CALIBRATION OF THE 
THEMATIC MAPPER AND MOLTISPECTRA.L SCANNER SYSTEM" 


Contract Number NASS- 27382 
For the Period: 1 November 1983 - 1 February 1984 


NASA/Goddard Space Flight Center 
Greenbelt MD 20771 


James M. Palmer, Co-Investigator 
Philip N. Slater, Principal Investigator 


Optical Sciences Center 
University of Arizona 
Tucson, Arizona 85721 






* — V 





Introduction 


§ 

% 

V- ■ 

I 

i 

I 

1 

v; 

$ 

I 

; 

I 

? 

: 



1 




1 




■i 




This is the fifth quarterly report on Contract NAS5-27382 entitled, 
"Spectroradioraetric Calibration of the Thematic Mapper and the 
Multispectral Scanner System." In this report, we provide the results of 
an analysis that relates TM saturation level to ground reflectance, 
calendar date, latitude, and atmospheric condition. We also include a 
revised version of the preprint included with the last quarterly report 
for publication in the IEEE Transactions on Geoscience and Remote 
Sensing. The paper is entitled 'Tn~flight Absolute Radiometric 
Calibration of the Thematic Mapper." 
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Sensor Radiance for a Midlatitude Atmospheric Model 

f 

To account for all orders of multiple scattering, the Herman j 

radiative transfer code iteratively traces solar radiance through a 
discrete number of scattering angles and many atmospheric layers. 

i 

Because of the large amount of computing time required, this has proved 
to be rather costly. Recently, however, the code has been transfered 
from the University computer system to an IBM personal computer. It can 
now be run continuously on this dedicated system. Although our system 
is equipped with 512K bytes of memory, only 192K bytes are required to 
run the code. A portion of the additional memory is configured as a 
"superdrive” where programs, such as the Fortran 77 compiler, are 
stored. Thus, when the source code needs to be recompiled, the code is 
read from RAM memory, rather than from floppy disks. Compilation time 
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is reduced by a factor of four, to about 10 minutes. Currently we are 
awaiting the arrival of Digital’s Fortran compiler (just released to 
dealers this week). This software will support complex arithmetic and 
also the 8087 math processor chip. It is believed that use of the 8087 
will reduce the run time of the Herman code by an order of magnitude. 
This would allow an average calculation to be done in about 15 minutes. 

The first data generated on the IBM are presented here. We wished 

to use the code to predict what ground reflectance would saturate the 

TM sensor output given a standard midlatitude atmospheric model. This 

model atmosphere is somewhat difficult to characterize due to the large 

variations in aerosol and water vapor levels. Although no single optical 

depth component can characterize aerosol scattering, measurements by 

Elterman 1 have allowed a rough correlation to be made between 

atmospheric visibility and the Mie optical depth, Tj^ e . At a wavelength 

of 0.55 um and with extremely clear atmospheric conditions a typical 

value of is 0.2. To bracket this aerosol level, values at 0.55 

um of 0.1 and 0.3 were also used as inputs. These three Tj^ e values can 

be correlated with atmospheric conditions corresponding to 20 km, 40 km, 

and greater that 100 km visibility. When using such a correlation, one 

must keep in mind that it is an extrapolation based upon data taken 

between what was defined as 2 km and 13 km visibilities. For any of 

these atmospheres a range of spectral values was computed using 

-<v-2) 

the relationship tHie ” bl . This spectral dependence was found by 

— { V+l ) 

assuming a Junge radial size distribution, with dN/dr => cr and 

v = 2.5. 




The ozone optical depth components were computed by using the U. S. 
Standard Atmosphere (USSA) Midiatitude Ozone Model value of 347 matm-cra 
for columnar ozone. 1 Using the spectral absorption coefficients of 
Vigroux, 5 t oz was computed for the wavelengths of interest. Absorption 
due to water vapor also had to be accounted for in band 4. The 
LOWTRAN V code was used to compute the transmittance of several lines 
within this band for the USSA water vapor content of 1.17 gm/cm l . These 
transmittances were averaged to compute a model x Hj0 value. Rayleigh 
scattering is accounted for by computing r Ra y using an atmospheric 
pressure of 1013 mb. Finally, the exo-atmospheric solar irradiance was 
accounted for by assuming the published values of Neckel and Labs,* 
given at an earth-sun distance of 1 AU. The values of these parameters 
are summarized in the tables to follow, as well as the code outputs. 

The output of the Herman radiative transfer code gives the 
upwelling radiance that is leaving the atmosphere and incident at the 
entrance pupils of the Landsat 4 and 5 sensors. This radiance is 
compared to the radiance level at which saturation occurs. The latter 
value is known from the prelaunch calibration in which a 122-cm 
integrating sphere was used to determine the gains, offsets, and dynamic 
ranges for each channel and each band. 5 Later, reference was made to 
the on-board internal calibration system to monitor po3tlaunch changes. 
Knowing the gain (G) and offset (0) coefficients, the maximum allowed 
radiance value can be computed by setting the digital counts (C) to 255 
in the calibration equation C = GL+0. This is the radiance which, when 
incident on the entrance pupil, would exceed 255 output digital counts 



following the use of the Radiometric Look-Up Table (RLUT). The 
saturation radiance values, as used in the 15 JAN 84 update to the TIPS 
data processing system are: 

HARD 1 2 3 4 5 7 

I. sa t 15.21 29.68 20.43 20.62 2.719 1.438 

(®W/ca* pa sr) 

Figures 1 through 4 plot the Herman code output as a function of 
solar zenith angle (from the data tabulated in Table 1). To utilize 
these figures, solar zenith angles at the time of Landsat overpass must 
be known. This is approximated in Figure 5, which gives these angles 
verses time of year and latitude.* The figure captions summarize the 
conditions required to saturate the various bands. Figure 6 shows that 
the ground reflectance corresponding to 255 digital output counts is 
relatively independent of aerosol loading. The last figure, Fig. 7, plots 
the radiance at the sensor for a ground reflectance of P ■ 0.75 and 
solar zenith angle of 45°. 
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Band 1 


X 

- 0.485 

T Mie 

- 0.2130 

Eo 

- 190.81 

T Ray 

- 0.1622 

L sac 

- 15.21 

T o* 

- 0.0100 



Fig. 1. Radiance at the sensor as a function of solar zenith angle for 
a visibility of 40 km. This shows, for example, that P“0.25 
does not give rise to saturation but p»0.50 is only unsaturated 
from mid-September to mid-March at 55 °N, from roughly mid- 
November to early February at 30 *N, and remains saturated 
through the whole year for 15 # N and 0 # N latitudes. 
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Band 2 
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Eq m 100.34 
^*tiat " 29.68 


“ 0.1965 
' r Ray “ 0.0839 
r oa " 0.0429 



Fig. 2. Band 2 remains unsaturated throughout the year for ground 
reflectances less than P « 0.50. The p =* 0.75 curve crosses 
the saturation line at a solar zenith angle of about 38°. Thus 
the sensor will output less than 255 digital counts from mid- 
July to May at 55°N latitude, from September to April at 30°N, 
from late September to March at 15°N, and from October to 
March at 0°N. 
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Fig. 3. In band 3 Che sensor remains unsaturated for a ground 
reflectance of p * 0.50 or less. At p=»0.75 the sensor is 
unsaturated from about early September to April at 55^ 
latitude, from mid-October through late February at 3t °N, and 
is saturated for lower latitudes for the entire year. 
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Band 4 
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» 0.84 

T Hie 

- 0.1618 

E 0 

* 104.19 

T Ray 

- 0.0175 

L aat 

*• 20.62 

T os 

- 0.0024; t 01 q ■ 0.0700 



25 35 45 55 65 

Solar Zenith Angle 


Fig. 4. Band 4 remains unsaturated for p“0,75 or less. At a ground 
reflectance of p a 1.00 the sensor is unsaturated from late 
July to mid-May at 55°N, from September to early April at 30°N, 
from early October to March at 15°N, and from early November 
to February at 0°N. 
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Fig, 6, Values of ground reflectance required to saturate the TM 
sensor for three different aerosol loading conditions 
(corresponding to about 20 km, 40 km, and greater than 100 km 
visibility) and a solar zenith angle of 45". The T^ie values 
are given for a wavelength of 0.55 um, and vary slightly from 
band to band. 
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Fig. 7. Radiance L(oiW/cm l sr Wa) for the first 4 TM bands, given a 
ground reflectance of p =* 0.75 and solar zenith angle of 
9 Z =* 45°. Values of T(^^ e are for X » 0.55 pm and correspond to 

3 levels of aerosol content (visibility). ! 
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Table I. Output radiance L(saW/cra* nr ua) from Heraan code. 


T Mie 


p - 0.25 

i> ** 0.50 

p aat 

0.1065 

25 

14.82 

27.06 

0.26 


45 

11.42 

20.63 

0.35 


65 

6.89 

11.85 

0.67 

0.2130 

25 

14.60 

26.26 

0.26 


45 

11.33 

19.92 

0.36 


65 

6.90 

11.30 

0.70 

0.3195 

25 

14.36 

25.36 

0.27 


45 

10.99 

1 9.1 8 

0.38 


65 

6.62 

10.75 

0.74 


Band 2 


T Mie 

0.0982 


0.1965 


L 


«* 

p ■ 0.50 

p - 0.75 

p sat 

25 

24.02 

35.76 

0.62 

45 

18.27 

27.07 

0.82 

65 

10.24 

14.96 

— — ' ” 

25 

23.08 

35.17 

0.64 

45 

17.65 

25.91 

0.86 

65 

9,53 

14.35 



25 

22.65 

33.99 

0.66 

45 

65 

16.93 

9.20 

25.32 

0.87 


0.2947 
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Band 3 


T Mle 

0.0913 


0.1826 


0.2739 


L 


0* 

P » 0.50 

P 0.75 

p aat 

25 

20.84 

31.35 

0.49 

45 

15.97 

24.00 

0.64 

65 

9.15 

13.64 

— — — — 

25 

20.22 

30.39 

0.51 

45 

15.34 

23.02 

0.67 

65 

8.60 

12.74 

— — — ■ - 

25 

19.90 

29.92 

0.51 

45 

15.01 

22.51 

0.68 

65 

8.31 

12.25 

— 
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Baad 4 

T Mie 

0.0809 


0.1618 


L 


°2 

P ■ 0.75 

P “ 1.00 

°sa£ 

25 

19.24 

25.83 

0.80 

45 

14.52 

19.49 

— 

65 

7.89 

10.56 

— — 

25 

18.98 

25.57 

0.81 

45 

14.22 

19.15 

— 

65 

7.57 

10.14 

— — - 

25 

18.67 

25.22 

0.82 

45 

13.89 

18.74 

— 

65 

7.25 

9.73 

— 


j 


0.2428 
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ABSTRACT 

Ground spectral reflectance and atmospheric spectral optical depth 
measurements made at White Sands, New Mexico on January 3, 1983, were 
used with an atmospheric radiative transfer program to determine the 
spectral radiance at the entrance pupil of the Landsat-4 Thematic Mapper 
(TM). A comparison with the output digital counts of the TM, when imaging 
the measured ground area, provided an absolute calibration for five 
detectors in TM bands 2, 3, and 4. By reference to an adjacent, larger 
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uniform area, the calibration was extended to all 16 detectors in each of 
the three bands. Pre-flight calibration results agreed with these in- 
flight measurements to 6.6%, 2.4 Z, and 12.9 2 in bands 2, 3, and 4 
respectively. 

I. IHTEODOCTIoa 

The Thematic Mapper (TM) multispectral scanner system was placed 
into earth orbit on July 16, 1982, as part of NASA's Landsat-4 payload. 

The entire system had been calibrated in an absolute sense at Santa 
Barbara Research Center before launch (1, 2). During flight, an internal 
calibration system monitors the calibration of the focal plane (3). To 
determine temporal changes of the absolute radiometric calibration of the 
entire system in flight, we initiated a program at White Sands, New Mexico, 
on January 3, 1983, to make spectroradiometric measurements of the ground 
and the atmosphere simultaneously with TM image acquisitions over that 
area (4). By entering our measured values into an atmospheric radiative 
transfer program (5), we determined the radiance levels at the entrance 
pupil of the TM in four of the TM spectral bands (fall width at half 
maximum: band 1, 0.45-0.52 urn; band 2, 0.53-0.61 um; band 3, 0.62-0.69 urn; 

and band 4, 0.78-0.91 urn). These levels were compared to the output 
digital counts from the detectors that sampled the radiometrically 
measured ground area, thus providing an absolute radiometric calibration 
of the entire TM system utilizing those detectors. Than, by reference to 
an adjacent, larger uniform area, the calibration was extended to all 16 
detectors in each of the three bands. 
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II. DKTKRMIHATIOH OP SPHCTEAL RikDIMCE AT THB TM 
FROM SITE MEASUREMENTS 

On January 3, 1983, an 80-mm layer of two-day-old snow covered the 
fl.at gypsum surface at the White Sands Missile Range. The reflectance of 
the snow was measured by reference to a 1.2 * 1.2-m barium sulfate panel 
using a Barnes Modular Multispectral 8-Channel Radiometer (Model 12-1000) 
(6), which collected radiant flux simultaneously in all the TM spectral 
bands over a total field angle of 15°. The Instrument was mounted on a 
rotatable boom 2.5 m above the snow to allow an average radiance value to 
be determined for an area of about 0.5 x 5 m. The measurements were made 
at 1708 GMT, coinciding with the overpass of the TM. The solar zenith 
angle was 62.8°. With the radiance of the barium sulfate panel at 62.8° as 
a reference, the reflectance of the snow in TM bands 1 to 4 was found to 
be 0.769, 0.761, 0.756, and 0.732 respectively, with an rms uncertainty of 
10.02 in all cases. 

After the TM overpass, hand-held Exotech (6) radiometer measurements 
were made at the pixel centers of a 4 x 4 pixel area on the snow (pixel 
size 30 x 30 m) to check for significant differences from the Barnes 
radiometer data. Variations between the Barnes and the Exotech data were 
small and attributed to variations introduced by the Exotech measurement 
technique rather than in the snow reflectance. The 4x4 pixel area was 
later found to be at an angle to the scan lines; as a result, the 
Individual detectors obtained different numbers of samples (see Table 6). 
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By reference to an access road and nearby frozen water surfaces we were 
able to accurately locate the known reflectance area on the TM Imagery. 

A solar radiometer (7) using nine 10-nm-wide spectral bands in the 
visible and near ir was used to determine the total spectral optical 
depths Tg XC at these nine wavelengths. The measured barometric pressure 
of 889.5 mbar allowed the Rayleigh spectral optical depth T' Ray to be 
determined. From these data the Mie and ozone spectral optical depths (8) 
were determined. 

In bands TM 1, 2, and 3, the component of x£ xt due to molecular 
absorption, x^ ba , is due entirely to ozone. In band 4, water vapor and C0 2 
predominate. Their effects have been included in Table 1, which lists the 
values for the various atmospheric constituents in bands 1 through 4. 

Because of the limitations of the instruments used on January 3, 
x'abs l n band 4 due to water vapor and CO * could not be measured. 
Instead, the LOWTRAN code (9) was used to compute the water vapor and C0 2 
transmittances across band 4 at 5 cm" 1 intervals. The water vapor 
transmittance was scaled, to account for the 44 % relative humidity 
measured at White Sands on January 3, then averaged to find t 'h 2 o+co 2 ' 
The effect of an error in this estimate can be judged by noting that the 
inclusion of water vapor and COj lowered the predicted radiance level at 
the TM by only 2%. 



TABLE 1. ATMOSPHERIC VALUES, 


Band 


T f 

T ezt 

T 'Mi« 

T *Bay 

T ’abB 

1 

0.485 

0.291 

0.148 

0.142 

0.001 

2 

0.57 

0.213 

0.138 

0.074 

0.006 

3 

0.66 

0.172 

0.128 

0.041 

0.003 

4 

0.84 

0.134 

0.110 

0.015 

0.009 


The data in Table 1 and Che reflectance values quoted above were 
used as inputs to the Herman and Browning (S) radiative transfer code. In 
using the code, the atmosphere is divided into a sufficient number of 
plane-parallel layers such that changes in radiance within each layer are 
due only to single scattering processes. 'the Gaus3-Seidel iterative 
technique is used to solve the equation of radiative transfer. Upon 
convergence, all multiple scattering effects have been taken into account. 
A Junge radial size distribution was assumed for the aerosols. A v value 
of 2.5 was assumed in the equation for radial distribution: 

SE . c,- (wl > (D 

dr 

where r is the aerosol radius. Values of 5.02 urn, 0.02 pm, and 0.04 pm 
were used for the maximum and minimum radii and incremental step size, 
respectively, for the aerosols. The aerosols were given a refractive 
index of 1.54-0.01 i and assumed to have a vertical distribution as 
measured by Elterman (10). 
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The following quantities were calculated by the code. Their values 
are listed in Table 2. 

Ed Dir : The downward direct solar irradiance at the ground is 

’ E«cos9 z exp(-Tg Xt sec0 z ) where E e is the exoatmospheric solar 

irradiance and 0 Z is the solar zenith angle 
Ep Dif s The downward diffuse solar irradiance at the ground 

Eu'dIt 2 The upward direct radiance at the TM due to Ed Dir ant * E D Dif ■ 

Ly|p: The upward path radiance at the TM is 

l T “ ( E D,Dir + E D,Dif) e3t P(“ T ext sec5 °) 

LpS The total radiance at the TM at a 5° nadir angle. 


TABLE 2. mSADIAHCES AMD RAIJIAHCES (NORMALIZED TO UHITT 
SOLAR EXOATMOSPHERIC IRRADIAHCE). 


BAUD 

SOLAR 

ZENITH 

AHGLE 

E D,Dir 

E. 

E D,Dif 

E. 

L 0,Bir 

8. 

L U,P 

8. 

L<j) 

17 

1 

1, 

i| 

tj 

t 

1 

55* 

0.345 

0.185 

0.097 

0.033 

0.130 



65* 

0.212 

0.152 

0.067 

0.025 

0.092 

: H 

2 

55* 

0.392 

0.145 

0.105 

0.024 

0.129 

* 


65* 

0.253 

0.122 

0.073 

0.019 

0.092 

■ : 

3 

55* 

0.425 

0.122 

0.111 

0.019 

0.130 

f 


65° 

0.282 

0.104 

0.078 

0.015 

0.093 

■ 

4 

55* 

0.454 

0.095 

0.112 

0.014 

0.126 

* 

j 

\ 


65* 

0.308 

0.082 

0.079 

0.011 

0.090 

< 

i 


Values for the exoatmospheric irradiances within the TM passbands, E 0 , 
were found by evaluating 
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K. 



( 2 ) 


' 1 * 

where d *» 0.983 astronomical units for January 3, X, and X a are the 
wavelength limits of the TM passbands as determined by using the 
equivalent passband technique of Palmer and Tomasko (11), and E 0 x Id the 
exoatmospheric solar spectral lrradiance data of Neckel and Labs (12). 
The integration was carried out using 2 to 10-nm intervals as dictated by 
the resolution of the Neckel and Labs data. The required values of L^/Eo 
for the center wavelength of each band and the solar zenith angle of 62.8° 
were found by interpolating between the 55“ and 65° radiative transfer 
code data in Table 2. These, when multiplied by their corresponding E 0 
values, gave the radiances in mWcnf^sr -1 in the TM passbands listed in 
Table 3. (An error of less than IZ is introduced using the value at the 
center wavelength of each band instead of evaluating 


E, 



oX‘*(X)d 



>X, 

where i(A) is the output as a function of wavelength /:rora the radiative 
transfer code, as given by L^/E, in Table 2.) 



TABLE 3. EX0ATMQSPEKRXC IRBADXA13CSS AM) TEE RADIANCES AT THE TH IS 
THE TM PASSBAHDS AS DBTEMIHID FEQM SITE HEASURBHBBTS AT WHITE 
SANDS. 



EQUIVALENT WAVELENGTH 

TM BAUD WIDTH LIMITS 

E e In 

L*j itt 

BAND 

in um 

IB usa 

mWcm~* 

bWc®' 

I 

0.0715 

0.4503—0*5218 

14.4 

1.45 

2 

0.0887 

0.5269-0.6156 

16.6 

1.66 

3 

0.0771 

0.6213-0.6984 

12.3 

1.25 

4 

0.1349 

0.7719-0.9068 

14.7 

1.44 


m. DETEEI3IHATIOS OF AVERAGE IH-BAHD RADIAHCE AT THE TM 
FROM IMAGE DIGITAL COOHTS AMD SYSTEM CALIBRATION DATA. 

We next describe how we determined corresponding sets of L-y values 
from pre-flight and in-flight internal calibration data. By identifying 
our site on the raw image data, we determined which detectors scanned 
the area and in what order, and how many samples each collected. The 
site was scanned from north to south by detector numbers 3, 2, 1, 16, and 
15, in that order. They collected 1, 3, 5, 4, and 2 samples respectively. 
We found the average digital count for each detector, using offset and 
gain values reported by Barker et al. (1), then calculated the average 
spectral radiance, L T , from L T * (average count-offset)/gain. (Throughout 
this paper, spectral radiances are referred to the entrance pupil). 
These spectral radiances were multiplied by the number of samples for 
each detector. The resultant products were added and then divided by 
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the total number of samples, 15. Titus we derived- a value for the 
average spectral radiance of our site as measured by the TM, 
proportionally weighted according to the number of samples per detector. 
These values are listed in Table 4. 


TABLE 4. AVERAGE COUNTS AMD SPECTRAL RADIAHCES (in raNc® - 2 er“* 
UB»“‘) FOR FIVE DETECTORS IE TM BABBS 2, 3 AMD 4 AS CALCULATED FROM 
PUS-PLIGHT CALIBRATION DATA (1). 



BAND 

2 

BAND 3 


BAND 4 


DETECTOR 

AVERAGE 

COUST 

average 
SPECTRAL 
KADIAK CB 

AVERAGE 

COUiJT 

AVERAGE 

SPECTRAL 

RADIANCE 

AVERAGE 

COUNT 

AVERAGE 

SPECTRAL 

RADIANCE 

3 

139.0 

17.04 

165.0 

15.40 

133.0 

11.89 

2 

143.7 

17.42 

169.7 

15.86 

132.0 

12.02 

1 . 

146.2 

17.52 

171.6 

15.73 

134.2 

12.00 

16 

141.5 

17.43 

167.8 

15.85 

132.0 

12.01 

15 

147.5 

17.65 

172.5 

15.88 

131.5 

12.05 

WEIGHTED 

AVERAGE 


17.46 


15.78 


12.01 


Then the average spectral radiance per band was multiplied by the 
equivalent bandwidth to provide the weighted average radiance in ' ach 
band. (These values are listed in column B of Table 5.) 

Barker et al. (2) have reported that the TM internal calibrator, 
used in flight, indicates that the response of TM bands 1 through 4 has 
slowly decreased with time during the period July-Deceraber 1982. We 
have used the Barker et al. data to change the pre-launch calibration 
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data. The changed values are Hated in column C of Table 5. With 
respect to this decrease in response, we note, first, that the internal 
calibrator compares the response of only the TM filters, detectors, and 
electronics to seven different irradiance levels; it does not measure any 
change in transmittance of the image-forming system. Second, it is 
possible that the decrease in response is wholly or partly due to a 
change in the output of the Internal calibrator. 

Table 5 summarizes our results baaed on measurements at White 
Sands and compares them to ■ pre-flight calibration and in-flight 
calibration data from the TM internal calibrator. 


f t$p****T 
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TABLE 5. COMPARISON OP RADIANCES (raWcia” •sr”') IN TM BAUDS 2, 3 AND 

4. 



COLUMN 

COLUMN 

COLUMN 



BAND 

A 

B 

C 

<T?>* 

<^r>* 

2 

1.66 

1.55 

1.50 

6.6 

9.8 

3 

1.25 

1.22 

1.16 

2.4 

7.2 

4 

1.44 

1.62 

1.57 

-12.9 

-9.6 


COLUMN A GIVES THE RADIANCE LEVELS AT THE TM AS DERIVED FROM 
GROUND Mm ATMOSPHERIC MEASUREMENTS AT WHITE SANDS ON JANUARY 3, 
1983, AMD THE USE OP AN ATMOSPHERIC RADIATIVE TRANSFER PROGRAM (SEE 
TABLE 3). 

COLUMN B GIVES THE WEIGHTED AVERAGE RADIANCE IH EACH PASS BAND 
(EiWcra~*8r~ ') AS DETERMINED FROM TA IMAGE DATA OF OUR WHITE SANDS 
SITE IH CONJUNCTION WITH PRE-FLIGHT CALIBRATION DATA. 

COLUMN C GIVES THE VALUES IN COLUMN B AS MODIFIED BY THE 
CHANGE m RESPONSE SUGGESTED BY THE INTERNAL CALIBRATOR DATA OF 
DECEMBER 8, 1982. 


Because of the small number of samples collected over the test 
site and the small variation in the digital counts recorded by each 
detector in each band, the most meaningful way to summarize the results 
of the calibration is to list the individual counts for each detector and 
the radiance in each band. This has been done in Table 6. Note that 
whereas there are significant detector-to-detector variations within a 
band, the variations for a given detector are consistent with the 

V 

expected uncertainty in the output of the analog-to-digital converter 
for a constant input analog signal. 
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TABLE 6. DIGITAL COUNTS FOE EACH DETECTOR THAT SAMPLED THE GROOM® 
SITE ABB THE RA2SIAECR AT THE TM DBTKBUBJiSD FROM THE HBASOOBMSHTS 
AT WHITE SAHBS. 








RADIANCE 

AT TM, 

DETECTOR 

BOMBER: 

3 

2 

1 

16 

15 

(nUcFs~ 2 or - ‘) 

BAND 2 

139 

143 

147 

141 

147 

1.66 

- 


144 

146 

141 

148 




143 

145 

141 






147 

143 






146 




BAND 3 

165 

170 

172 

167 

172 

1.25 



169 

172 

167 

172 




170 

171 

168 






171 

169 






172 




BAND 4 

133 

132 

134 

133 

130 

1.44 



131 

134 

131 

133 




133 

134 

131 






134 

133 




135 


IV. EXTENSION OP CALIBRATION TO OTHER DETECTORS 


To extend the calibration of the five detectors in each band to all 


16 detectors in each band, a 10 « 16 pixel test site of uniform 
reflectance was selected from an analysis of the TM imagery of White 





Sands. This site was 500 m southwest of the 4x4 pixel test site 
described earlier. Its uniformity was such that none of the 16 rows of 
10 pixels exhibited a variance of greater than 1.5 digital counts, the 
average variance being 0.9 digital count. The average number of digital 
counts for each of the 16 detectors in each band in this 10 x 16 pixel 
area is listed in Table 7. To verify, a posteriori, that the large site 
had the same reflectance as the 4x4 pixel site, the digital counts 
were compared for detectors 3, 2, 1, 16, and 15 in the three bands. For 
the 15 detectors that sampled both sites (five each in three bands) the 
digital count difference between the two sites was less than 1 for nine 
of the detectors. 

It is interesting to note that there was no evidence of a level 
shift, as described by Kieffer et al. (14), in this part of the scene. The 
10 * 16 pixel area was selected on the basis of its close statistical 
match to the 4x4 pixel area and was later found to have been scanned 
by detectors 1 through 8 in one direction and by detectors 9 through 16 
In the reverse direction. However, in other parts of the scene, certain 
detectors did exhibit a level shift of up to 3 counts. 

The average digital count values in Table 7 can be compared to our 
computed TM entrance pupil radiance values given in Column A of Table 5 
to provide a single point calibration for all detectors in TM bands 2, 3, 


and 4. 


- • * -'m • 
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TABLE 7. AVERAGE DIGITAL COUNTS FOE THUS 10 x 16 PIXEL TEST SITE IN 
TEE THESE TM BAUDS. 


DETECTOR 

BAND 2 

BAND 3 

BAND 4 

1 

145.4 

171.5 

133.8 

2 

142.7 

167.8 

131.0 

3 

141.0 

167.4 

133.0 

4 

143.0 

166.8 

131.4 

5 

145.2 

168.4 

131.2 

6 

144.0 

168.7 

132.4 

7 

142.3 

167.3 

131.7 

8 

142.5 

167.5 

132.0 

9 

143.4 

168.3 

132.6 

10 

142.1 

167.5 

137.2 

11 

141.4 

166.3 

130.9 

12 

141.7 

169.8 

132.7 

13 

145.2 

168.6 

134.2 

14 

143.8 

168.9 

132.3 

15 

146.2 

170.2 

130.7 

16 

141.0 

167.1 

132.1 


V. UNCERTAINTIES IN THE CALIBRATION PROCEDURE 
There are four sources of uncertainty in the calibration procedure. 
These have been discussed elsewhere (4) and will be only briefly referred 
to here. 

The main uncertainty is in our knowledge of the composition, size 
distribution, and vertical distribution of atmospheric aerosols. Of these, 
the value assumed for the aerosol complex index is the most critical. On 
the basis of a survey of aerosol conditions at White Sands by Jennings et 
al. (15), we estimated that an error of ±3.5% is introduced if extreme 
rather than the assumed average conditions existed on the morning of the 
calibration measurements. Departures of aerosol size and vertical 
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distributions from the assumed average can give rise to a ±2% change in 
radiance at the TM. We thus have a root-sum-of-the-squares (rss) 
uncertainty of *4 % due to an imprecise knowledge of atmospheric 
conditions. 

An uncertainty of ±2 % is attributed to the combination of three 
reflectance-related factors: nonuniformities in the ground reflectance, 
uncertainty in the knowledge of the absolute reflectance of the reference 
panel, and the departure of the ground and reference surfaces from being 
lambertian reflectors. 

An uncertainty of ±1 % is due to the uncertainty in the 
exoatmospheric spectral irradiance data of Neckel and Labs (11). 

An uncertainty of 10.7% is due to a variation of H digital count in 
an output of 150 counts from the TM in the imaging of a constant 
radiance level. 

If we assume an error of ±ZZ due to the neglect of polarization 
effects, the overall rss uncertainty in the calibration procedure is ±5%. 

Henceforth, to reduce the uncertainty t:o less than ±3%, we plan to 
take the following steps: 

1. Collect aircraft 3pectropolarimeter data, calibrated in an absolute 
sense, to check the calibration procedure and determine the 
magnitude of any possible systematic errors. 

2. Collect aircraft imagery of the ground reference areas to 
accurately map spatial variations in ground reflectance. 
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3. Make ground measurements simultaneously with TM overpasses at 
sites having different reflectances* 

4. Use a version of the Herman and Browning atmospheric radiative 
transfer program (5) that accounts for polarization and the 
nonlambertian characteristics of the ground surfaces. 

5. Attempt to more precisely estimate the aerosol composition and 
distributions by: 

a. Making ground level diffuse and direct irradiance 
measurements. 

b. Applying inversion techniques. 

c. Making more detailed atmospheric measurements, e.g. of the 
Babinet and Arago points. 


VI. CONCLUSIONS 

The absolute calibration of five detectors in TM bands 2, 3, and 4, as 
determined by measurements at White Sands on January 3, 1383, is given in 
Table 6. Pre-flight calibration results agree with these in-flight 
measurements to 6.6%, 2.4% and 12.9% in bands 2, 3, and 4 respectively. 
Table 7 shows the absolute calibration extended to all the detectors in 
bands 2, 3, and 4. The estimated uncertainty in these results is ±5%; the 
estimated uncertainty in the pre-flight calibration is no better than ±6% 
(1. 16). 
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Band 1 saturated over the snowfield at White Sands. Pre-flight data 
(1) indicate that a saturation level of 255 counts corresponds to a 
radiance at the sensor of 1.14 mWcra~ 1 ar“ l in TM band 1. We estimate that 
the snowfield provided a radiance level of 1.45 mWcm“ , 3r“ 1 at the sensor. 

We are presently fabricating field equipment to provide more detailed 
and accurate measurements of the surface and atmosphere at White Sands. 
Our goal is to reduce the uncertainty in sensor absolute calibration to 
less than ±3%. 

We plan to continue the work described here to include the in-flight 
absolute radiometric calibration of the Landsat-5 Multispectral Scanner 
System and TM and the Systeme Probatoire d'Observation de la Terre, Haute 
Resolution Visible (SPOT/HRV) systems. 
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